arXiv:1509.03330v2 [hep-ph] 5 Nov 2015 


TTP15-034 


Integration by 
A criterion 


parts identities in integer numbers of dimensions, 
for decoupling systems of differential equations 


Lorenzo Tancredi*’’ 0 


^ Institute for Theoretical Particle Physics, KIT, Karlsruhe, Germany 


Abstract 

Integration by parts identities (IBPs) can be used to express large numbers of apparently different 
d-dimensional Feynman Integrals in terms of a small subset of so-called master integrals (Mis). Using 
the IBPs one can moreover show that the Mis fulfil linear systems of coupled differential equations in the 
external invariants. With the increase in number of loops and external legs, one is left in general with 
an increasing number of Mis and consequently also with an increasing number of coupled differential 
equations, which can turn out to be very difficult to solve. In this paper we show how studying the 
IBPs in fixed integer numbers of dimension d = n with n € N one can extract the information useful 
to determine a new basis of Mis, whose differential equations decouple as d —>■ n and can therefore be 
more easily solved as Laurent expansion in (d — n). 
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1 Introduction 


Dimensionally regularised m Feynman integrals fulfil different identities, among which the most notable 
ones are the so-called integration by parts identities (IBPs) [11[5]. Given a family of Feynman integrals, 
the IBPs can be used to write down a large system of linear equations with rational coefficients, which 
contain the Feynman integrals of that family as unknowns. By solving algebraically the system a large 
number of apparently different Feynman integrals can be expressed in terms of a much smaller basis of 
independent integrals dubbed master integrals (Mis). In realistic applications the number of such equations 
can grow very fast, requiring the use of computer algebra in order to handle the complexity of the resulting 
expressions. There are different public and private implementations which allow to perform the reduction 
to Mis in a completely automated way [Ma] based on the so-called Laporta algorithm [lOllTl] . 

The IBPs can be used to prove that dimensionally regularised Feynman integrals fulfil linear systems 
of first order differential equations in the external invariants [IMS]- A thorough review of the method 
can be found in m- Considering a Feynman graph which is reduced to N independent Mis, by direct use 
of the IBPs one can derive a system of N coupled linear first order differential equations for the latter, 
which can be rephrased as an A^-th order differential equation for any of the Mis. Supplemented with N 
independent boundary conditions, the system of differential equations contains all information needed for 
numerical or analytical calculations of the Mis. Indeed, in the general case, the analytical solution of an 
iV-th order differential equation is a very non-trivial mathematical problem. 

It has been observed that, in many cases of practical interest, a substantial simplification of the problem 
occurs when studying the behaviour of the system of differential equations as the space-time dimension 
parameter d approaches 4, which is also the physically relevant case. Usually we are indeed not interested 
in an exact solution for the Mis as functions of d, but instead in the coefficients of their Laurent expansion 
for d ss 4. In |18lll9j . and in many subsequent applications of the differential equation method, it was shown 
that it is often possible to choose a basis of Mis such that the differential equations take a simpler triangular 
form in the limit d —^ 4. If this is possible, the problem of integrating the system of differential equations 
simplifies substantially, reducing de facto, at every order in (d — 4), to N subsequent integrations by 
quadrature. Experience showed that, whenever such a form is achievable, the differential equations can be 
integrated in terms of a particular class of special functions, the multiple polylogarithms (MPLs) [T81I20II21) . 
The latter have been studied extensively by both mathematicians and physicists and routines for their 
fast and precise numerical evaluation are available since some time Disclosing their algebraic 

properties allowed furthermore the development of very powerful tools for the analytical manipulations of 
these functions [25H?7] . 

More recently it has been shown [281130] that in many of these cases a basis of Mis can be found, such that 
the system of differential equations takes a particularly simple form, commonly referred to as canonical 
form. The system is said to be in canonical form if the regularisation parameter d can be completely 
factorised from the kinematics, appearing as an explicit (d — 4) factor in front of the matrix of the system. 
In addition, the coefficients of the matrix must be total differentials of logarithms of functions of the external 
invariants (i.e. they are said to be in d-log form). A canonical basis is particularly convenient as it allows 
a straightforward integration as series expansion in (d — 4) and, due to the d-log form of the coefficients, 
it integrates directly to MPLs of uniform transcendental weight. Criteria for the construction of candidate 
canonical integrals have been presented in [29] and developed in detail, for example, in [31] . In the special 
cases in which the differential equations depend only linearly on the dimensions d, the Magnus algorithm 
can be used to perform a rotation of the system to a canonical form [32] . For a recent application of the 
algorithm see [33]. A completely different approach based on Moser algorithm [34] has been developed for 
the univariate case in [35] and discussed also independently in [36] . Another interesting approach is based 
on the properties of higher order differential equations fulfilled by the individual master integrals [37] . In 
spite of all this impressive progress, a fully automated algorithm, working also in the multivariate case, 
is still missing. It has nevertheless become clear that, if one can find a basis of Mis whose differential 
equations become triangular as d —>■ 4, it is often (but not always) relatively easy to bring it in canonical 
form by removing the undesired terms in the differential equations [38] . 

Indeed, different examples are known where neither finding a canonical basis nor a triangular one is 
possible. In all these cases the master integrals cannot be integrated in terms of simple multiple poly- 
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logarithms only [5^42) . It becomes therefore a very interesting problem that of finding a set of criteria 
to determine whether, given a Feynman graph, there exists a basis of Mis whose system of differential 
equations becomes triangular in the limit for d —>■ 4. This would provide a way to easily classify all possible 
diagrams where this is not possible and that would therefore be expected to evaluate to more complicated 
classes of function. Aim of this paper is to show that a large amount of information about the possibility 
of achieving such decoupling can be extracted by studying the IBPs for fixed even numbers of dimensions, 
i.e. in the limit d —^ 2 n, where n is any integer number. Tarasov-Lee shifting identities [431144] allow, in 
fact, to directly relate the structure of the differential equations in any even number of dimensions with 
the physical case d = 4. Indeed, since Feynman integrals are usually divergent at d = 2 n, this limiting 
procedure must be carried out as a Laurent series in (d — 2 n). 

The rest of the paper is organised as follows. In Section [5] we review the use of integration by parts 
identities for the reduction to master integrals and we summarise the main results of the differential equa¬ 
tions method. In Section [3] we outline the central idea of the paper. We show in particular what kind 
of information can be extracted by solving the IBPs as Laurent series for d —>■ 2 n and how to use it to 
simplify the system of differential equations. In Section |4] we then apply these ideas explicitly to many 
different examples of increasing complexity. Some comments on the method are given in Section 14.71 where 
we try as well to point out some relevant open issues. Finally we conclude in Section |SJ In Appendix lAl we 
compare our method with the Schouten identities introduced in m, while in Appendix |B] we show how to 
shift a system of differential equations of an even number of dimensions. 


2 Integration by parts identities and differential equations 


Let us consider an Z-loop scalar Feynman integral depending on P independent external momenta pi 


I{d; ai, ...,ar,bi, ...,ba) 


d% 

J bf (27r)2 Dl\..Dp 


( 2 . 1 ) 


where ki are the loop momenta, Di = — ml are the propagators and Si = kj ■ pi are irreducible scalar 

products. In view of the discussion below, we will usually denote the integrals of a topology as in (|2.1L 
keeping explicitly only the dependence on the dimensions d and on the powers of denominators and scalar 
products. In dimensional regularisation, for every integral of the form (12.111 there always exists a value of 
the space-time dimensions d, such that the integral is convergeni@. Necessary condition for the convergence 
of an integral is the integrand be zero at the boundaries. This condition can be mathematically rephrased 
as 



= 0 , 


( 2 . 2 ) 


1=1 


where the v!^ are any of the external or internal momenta = {fci,..., ...,pp} . Such an identity is 

called an integration by part identity or IBP. It is clear that in this way l{l + P) IBPs can be established 
for each integrand. Upon explicitly evaluating the derivatives and contracting with the momenta v^, new 
integrals belonging to the same topology (i.e. integrals with the same set of denominators) are generated. 
In particular, each IBP identity can relate integrals with (s — 1), s and (s -I- 1) powers of scalar products, 
and {t + r) or {t + r + 1) powers of propagators. Notice that, by contracting with new reducible 
scalar products can be generated, which could then simplify some of the denominators producing integrals 
belonging to any of the (t — 1) sub-topologies of the original graph. 

It has been shown [iniiiiiiii] that the system of IBPs, which appears in general to be over-constrained, 
can instead be solved allowing to express most of the integrals as linear combinations of a small subset of 
basic integrals, dubbed master integrals (Mis). Indeed, as for any algebraic basis, the choice is not unique 


2 All scaleless integrals in dimensional regularisation are zero for consistency. 
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and, by suitably changing the basis, one can substantially simplify the calculation of the integrals, as we 
will discuss later in this paper@ 

Let us consider now a Feynman graph (or a topology) characterised by a set of propagators Di and 
irreducible scalar products Si. All integrals will depend on the space-time dimensions d and on the external 
invariants Xij = Pi ■ pj, where Pi are as usual the external momenta. Let us assume, for definiteness, that 
by solving the system of IBPq all integrals for the given graph can be expressed in terms of a basis of N 
independent Mis Ii{d] Xij) with i = 1,N, which are of the same form of Eq. (12.11) . Differentiating with 
respect to any of the external invariants Xij amounts to differentiating with respect to linear combinations 
of the external momenta pf [14]. Therefore, by acting with these differential operators directly on the 
integrands of m, one produces linear combinations of integrals belonging to the same Feynman graph 
and to its sub-topologies. The latter can again be reduced to Mis, generating in this way a system of N 
linear first order differential equations with rational coefficients in any of the invariants Xij. Suppressing the 
dependence on the sub-topologies, which can be considered as a known inhomogeneous term in a bottom-up 
approach, the homogeneous part of the system can always be written as 


d 

dxij 


Xi (d, Xij ) 
(d, Xij ) 


eii{d^Xij') ... (^d^ Xij') 

C7vi(d, ... cpj]\j{d^ Xij) 


Ti (d, Xij ) 
Xtv (d, Xij ) 


(2.3) 


where the coefficients Cij{d;p^) are simple rational functions of the dimensions d and of the external invari¬ 
ants Xij. We can rewrite the system in matrix form as 


-^^I{d-,Xij) = A{d-,Xij)i{d]Xij), (2.4) 


where we introduced the vector of master integrals I{d]Xij) and the matrix of the coefficients A(d]Xij). 

The system (|2.3p is, in the general case, coupled and can therefore be rephrased as an 7V-th order 
differential equation for any of the Mis I,; (d; Xij). In most practical applications we are interested in 
determining the Mis as Laurent expansion for d « 4 

OO 

I.(d;x.,)= ^ x(“)(4;x.,)(d-4)“. (2.5) 

OL— — a 

By expanding both left- and right-hand side of ()2.3p one is left with a chained system of N differential 
equations where, at any order a, the previous orders can only appear as inhomogeneous terms. 

2.1 An optimal basis of master integrals 

It has been shown by Tarasov and Lee |43||44] that the value of a Feynman integral in d space-time 
dimensions can be directly related to that of the same integral in d — 2 or d -I- 2 space-time dimensions. This 
implies that, if all Mis of a given graph are known as Laurent expansion in any even number of dimensions, 
d = 2 n. 


OO 

li{d;xij)= ^ l|“^(2n;xij) (d-2n)“ , (2.6) 

a.— — b 

then the coefficients of their series expansions in d = 4, (4; x^) in (12.511 . can be obtained as linear 

combinations of the (2 n; x^). For more details see for example |39l41j and the discussion in ApoendixlBl 

different approach to reduction to Mis using hyperelliptic curves, and its equivalence to the IBPs in some explicit cases, 
was recently discussed in m- 

'^Note that, in order to have a complete reduction, one must also consider all possible symmetry relations among the 
integrals due to shifts of the loop momenta. 
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Indeed, changing the basis of Mis changes the form of the matrix A{d\Xij) in equation (12.41) . An 
interesting problem is therefore how to define an optimal basis of Mis in order to simplify as much as 
possible the system (lOl) . Since we are interested in computing the Mis as Laurent expansion in (d —4) (or, 
in general, in (d — 2 n)), an obvious simplification would occur if we could decouple some of the differential 
equations, at least in the considered limit. In particular, given a system of N coupled equations, one 
could think of classifying the complexity of the latter by determining the minimum number of differential 
equations that cannot be decoupled in the limit d —>■ 4 (or more generally d —>■ 2 n). 

At this point it is useful clarify more precisely what we mean with decoupling in this context. Let us 
consider a 2 x 2 coupled system of differential equation^ 

9 -* 

— I(d; x) = A(d; x) I(d; x ), (2-7) 

ox 

where I{d]x) = (/i(d; x),/ 2 (d; a;)) is the 2-vector of unknown functions, A{d,x) is a 2 x 2 matrix, d are 
the space-time dimensions and x is a variable the functions depend or@. Assume now for simplicity that 
the functions /i(d; x),/ 2 (d; x) are finite in the limit d —>■ 4 and that the matrix A(d; x) does not contain 
any explicit poles in l/(d — 4). Assume finally that, in the limit d —>■ 4, the matrix A(d;x) has non-zero 
non-diagonal entries, and therefore that the system is coupled in the limit d —> 4. Of course, by expanding 
the entries of A(d; x) in (d — 4) we can write our system as 

■^I(d;x) = A(°U4;x)X(d;x)-h(d-4)A^^h4;x)I(d;x)-hO((d-4f) . (2.8) 

ax ^ ^ 

It is now clear that if, by any means, we can find two independent solutions to the 2x2 system 

^ /(x) = A^°'> (4; x) /(x), (2.9) 

say (ni(x), U 2 (x)) and {wi{x),W 2 (x)), then we can define the new vector J(d]x) through the rotation 

J(d;x) = G(x)f(d;x), with G(x) = ^ W2(^) ) ’ 

such that the differential equations satisfied by J (d; x) assume the form 

^ J(d;x) = (d-4)G-i(x)A(i)(4;x)G(x)J(d;x) + G((d-4)2) , (2.11) 

i.e. they become trivial in the limit d —>■ 4. The matrix G(x) can be of course arbitrarily complicated, as it 
contains the solutions of a second order differential equation. In this case we would have of course achieved 
a decoupling, but at the price of having to solve a coupled system of differential equations, which is in the 
general case not possible. 

On the other hand, what we are really interested in is to determine whether a basis of Mis exists such 
that some of the non-diagonal terms of the matrix A(d; x) become zero in the limit d —>■ 4, and such that 
this basis can still be reached from our starting basis only through IBPs (i.e. without having to solve a 
coupled system of differential equations!). What this means in practice is that, if such a basis existed, then 
the rotation matrix G would assume a very simple form, namely it would contain only rational functions 
of the external invariants x^ (and of the dimensions d). This new basis would therefore fulfil a system of 
differential equations where some (or all) of the Mis decouple in the limit d —>■ 4, and still it would be a 
system of linear differential equations with rational coefficients only. In this respect we note that, for all 
known cases of Mis which can be integrated in terms of multiple polylogarithms, a change of basis in the 
sense described above can be found and the system of differential equations can be put in triangular form 
as d —>■ 4 

^ J(d;xy) = r(4;xy )y(d;xy ) + 0(d - 4), (2.12) 

^Again, we neglect the inhomogeneous terms everywhere. 

®In the case of Feynman integrals x represents a generic Mandelstam variable. 
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where T{4:;Xij) is a triangular matrix and does not depend on the dimensions d. From the point of view 
of the classification outlined above this corresponds to the easiest case, where all equations decouple in the 
limit d —>■ 4 and, effectively, the problem reduces to a series of independent integrations by quadrature. 
Finding a basis in this form is often a first step towards a canonical basis in the sense introduced in |29) . 
We want to stress here that of course all these considerations apply in the very same way for any integer 
number of dimensions d —>■ n (even or odd). 

Unfortunately a change of basis of this kind cannot always be found. Several cases are known where the 
system cannot be completely triangularised and instead at least two differential equations remain coupled. 
In these cases MPLs turn out not to be enough for describing the solution and the class of functions must 
be enlarged to include also elliptic generalisations of the latter [4^. It is unclear whether this will be the 
end of the story, since cases where three or more coupled equations survive are relatively easy to find, as we 
will show in the following. What still appears to be missing is a criterion to determine, given a Feynman 
graph, what is the minimum number of equations which cannot be decoupled. Together with simplifying 
as much as possible the problem at hand, this could also give a hint to which class of functions are required 
for describing the solution. 

3 Reading the IBPs in fixed numbers of dimensions 

In order to find a possible working criterion to determine the minimum number of coupled differential 
equations we should go back to think how the differential equations are derived. We saw that differentiating 
a master integral with respect to the external invariants produces new integrals belonging to the same 
Feynman graph. By using the IBPs one can then reduce these integrals to Mis, ending up with a system 
of differential equations. If we start with N master integrals we will obtain in general a coupled system of 
N linear differential equations. The fact that the N differential equations are coupled can be seen, in this 
respect, as due to the linear independence of the N master integrals in d dimensions. 

As we already discussed, for any physical application we are interested in computing Feynman integrals 
as Laurent series in {d — 4) or, more generally, in {d — 2 n) with n S N. Of course, different integrals 
have different degrees of divergence, i.e. their Laurent expansion starts at different orders in {d — 2 n). 
For any value of the dimensions, nevertheless, the maximal divergence can be computed in dimensional 
regularisation and depends only on the topology of the graph under consideration (i.e. on the number of 
loops, of external legs etc.). We can therefore imagine to first generate the IBPs in d dimensions and then 
expand them as a Laurent series in [d — 2 n), obtaining in this way a chained set of systems of IBPs, one 
for every order in (d — 2n). It is clear that, by construction, at every order in {d — 2n), the homogeneous 
part of each system will be identical, while the inhomogeneous part will contain the previous orders of the 
expansion (and the sub-topologies, that we will neglect throughout). If we limit ourselves to the first order 
of the expansion, i.e. the one corresponding to the highest pole in {d — 2n), the system of equations that 
we are left with is equivalent to the original system of IBPs where d is fixed to be d = 2 n, and corresponds 
to the sole homogeneous system. 

Now, it is very well known that upon fixing the number of space-time dimensions in the IBPs to an 
integer value it may happen that some of the equations degenerate and, in particular, that some of the 
integrals that used to be linearly independent for generic values of d, become linearly dependent from each 
other. From the point of view of the differential equations satisfied by the master integrals, if some of the 
integrals were to become linearly dependent in the limit d —^ 2 n, one would expect that those masters 
should not bring any new information in that limit and it should therefore be possible to decouple them 
from the system of differential equations as d —^ 2 n. Let us try to state this point more precisely. As 
exemplification we consider a topology that is reduced to 2 master integrals which we call Ii(d;x) and 
l 2 (d; x), where d are the dimensions and x is a generic Mandelstam variable. Neglecting the sub-topologies 
the system of differential equations that they satisfy can be written as 

{ ^Ii(d;x) = cii(d;x)Ii(d;x) -I- Ci 2 (d;x)l 2 (d;x) 

(3.1) 

^l 2 (d; x) = C 21 (d; x) Ii (d; x) -I- C 22 (d; x) I 2 (d; x) 
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where the Cij{d]x) are rational functions. Let us now follow the argument above and generate the IBPs 
fixing d = 2n. Let us assume that, by solving this simplified system, one of the two master integrals 
becomes linearly dependent from the other one and the new IBPs produce the relation 

X2{2 11 ] x) = b{x) Xi[2 11 ] x) ^ (3-2) 

where b{x) is a rational function of the Mandelstam variable^ Equation (13.21) implies that in d = 2 n 
one of the two master integrals becomes linearly dependent in the sense of the IBPs. According to the 
argument above we would therefore expect to be able to decouple the two differential equations in this 
limit. In order to see this it is useful to ask ourselves how such a relation can emerge from the original 
d-dimensional IBPs. Let us imagine that upon solving the IBPs for generic d, we can find a d-dimensional 
relation expressing a given integral of the graph under consideration, say K{d]x), as a linear combination 
of the two masters and such that 

K{d,x)= ^ {bi{d]x)Xi{d]x)+ b 2 {d]x)X 2 {d]x)) , (3.3) 

d — 2n 

with bi{x)/b 2 (x) = b(x) and n bi{d] x) = bi{x), for z = 1,2. It is clear that, if this is the case, the 

IBPs which would generate this identity for generic d, would instead generate ()3.2I1 once d is fixed to be 
d = 2n. These relations are precisely what we are looking for. To refer to the latter we will often use 
throughout the paper the notation 

bi{d] x)Xi{d] x) & 2 (d; x)X 2 {d] x) = 0{d — 2 n), (3.4) 


or equivalently 

bi{x)Xi{d]x) + b 2 {x)X 2 {d]x) = 0{d — 2n), (3.5) 

where it should be understood that, in general, this does not mean that the combination above is really of 
order 0{d — 2 n), but simply that it becomes zero upon setting d = 2n in the IBPs. Note that, of course, 
using the bi{x) instead of the bi{d]x) can only produce corrections of order 0{d — 2n) due to (13.3|) . We 
will see many examples of these relations in the sections below. 

Naively, the fact that only one integral is linearly independent for d = 2n would require that the integral 
itself should satisfy a first order differential equation as d —>■ 2 n. Finding a basis of master integrals for 
which this first order equation emerges is equivalent to finding a basis which decouples the system dSH. 
To this aim let us perform the following rotation of the master integral basis 


Ji{d]x) = bi{x)Xi{d]x) + b 2 {x)X 2 {d] x ), J 2 (d] x) = X 2 {d] x). 

The system dSlD under this rotation becomes 


(3.6) 


d 

— Ji (d; x) = ( cii (d; x) + 


b 2 {x)c 2 i{d] x) + b'i{x) 
bi{x) 


Ji{d]x) 


+ ( 6 i(a;)ci 2 (d; x) + 62 ( 2 ^) (c 22 (d; x) — cii(d; x)) + & 2 (x) — 


hjx) [b 2 ix)c 2 i{d]x) +b[{x)] 
bi{x) 


J2{d]x) 


■ J2{d]x) = 


,, . hix) , . 

C22{d]x) - 7—-^C2i(d; x) 

bi{x) 


J2{d]x) 


C 2 i(d; x) 
bi(x) 


Ji{d]x). 


(3.7) 


Equations (1X71) do not look particularly illuminating at first glance. We claim nevertheless that these 
equations are precisely what we were looking for. The basis Ji{d]x), J 2 (d]x) defined in (13.61) . in fact, has 
been chosen in order to exploit the linear dependence of the two master integrals in the limit d —> 2 n. In this 
limit the IBPs tell us that Ji{d] x) is by construction suppressed by a factor (d—2n) and therefore decouples 

^To be precise we should recall that, since the master integrals can be divergent, this relation cannot be seen, in general, 
as a real relation between the two masters. 
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from the problem. We expect therefore that the differential equation for the latter should decouple in this 
limit or, in other words, that 


(^i{x)ci 2 {d\x) + 62 ( 2 ;) {c 22 {d;x) - cii(d;x)) + 63 ( 2 :) 


b 2 {x) [b 2 {x)c 2 i{d;x)+b[{x)] \ ^ _ 2 

bi{x) ) 


(3.8) 

If this is true then upon expanding the system of differential equations as Laurent series in (d — 2 n) one 
can, at every oder, first solve the differential equation for Si{d]x) by quadrature, and then use this as an 
input for the second equation. A rigorous mathematical proof of equation (13.81) is outside the scope of this 
paper and we will limit ourselves to show explicitly how this works in practice with several examples of 
different complexity. 

The considerations above can be easily generalised to N master integrals Ii,..., Xjy. In this case one 
starts with a system of N coupled differential equations. By solving the IBPs for d = 2n one can then 
verify how many of the master integrals become linearly dependent in this limit. Assuming that N — M 
integrals remain independent, this means that M relations like (13.3|) can be found, say 


Ki{d;x) =-—— ( 6 ii(d;x)Ii(d;a;) + ... + bijv(d; x)Zjv(d; x)) 
a — 2 n 


KM(d;x)=-—— (bMi(d;x)Zi(d;x) + ... + bMN(d;x)Zjv(d;x)) , (3.9) 

a — 2n 

and the bij{d]x) are as always rational functions of the dimensions and of the Mandelstam variable^. As 
for the previous example we will often write these relations as 

bii{d] x)Zi{d- x) + ... + biN{d',x)Z]\f{d-x) = 0{d— 2n) 


bMi{d] x)Zi{d- x) + ... + bMN{d', x)ZN{d] x) = 0{d — 2n), (3.10) 

where once more we imply that these combinations become zero upon setting d = 2 n in the IBPs. As 
before we define bij{x) = limd ^2 n bij{d; x) and, following the same reasoning, we can then try to rotate 
the basis of master integrals to 


f7i(d; x) = bii{x)Zi{d; x) + ... + &iAr(x)lAr(d; x) 

JM{d;x) = bMi{x)Zi{d]x) + ... + bMNix)ZN{d;x) 

JM+iid;x) =ZM+iid;x) 

JN{d;x) =ZN{d;x). (3-11) 

Under the rotation (j3.11ll . we expect the M integrals f7i(d; a:), ..., JM{d', x) to decouple from the remaining 
independent integrals in the limit d —>■ 2n, as in (13.7L One must be cautious here on what is intended 
by decoupling. According to the arguments above, upon the change of basis (13.111) , we expect the system 
of differential equations to split into two blocks in the limit d —)• 2n, one M x M and the other {N — 
M) X [N — M). This would correspond, order by order in (d — 2n), to an M-th plus an (TV — M)-th 
order differential equation, unless for some other reason internally the two blocks of differential equations 
further decouple in this limit. On the other hand, for the Feynman graphs that we considered so far (see 
for example Sections 14.3l and r4.6|) . even a stronger claim can be made. In these cases the rotation (I3.11|) not 
only splits the system into two blocks, as described above, but it also produces an explicit (d — 2 n) in front 
of the whole M x M block originating from relations (13.101) . This explicit overall factor allows to effectively 

®For this to be true the relations must be linearly independent in the limit d ^ 2n. 


7 














reduce the problem to the solution of one single {N — M)-th differential equation, plus M integrations by 
quadrature. The reason for this behaviour is still partly unclear and deserves further study. 

Summarising, the discussion above brings us to the following conclusion. Given a topology with N 
master integrals which fulfil a set of N coupled differential equations in d space-time dimensions, the study 
of the IBPs in fixed numbers of dimensions, say d = 2n, provides a tool to determine how many master 
integrals can be decoupled from the differential equations as d —?> 2 n. Of course, the arguments given above 
are partly oversimplified and we have not provided here any rigorous mathematical proof. The structure 
of the differential equations can be in general very involved and, instead of embarking on complicated 
mathematical arguments, we prefer to show explicitly how this ideas can be simply applied to different 
cases of increasing complexity. In the next section we will start off by considering simple examples where, 
by fixing the number of dimensions to an even integer value, only one master integral remains linearly 
independent and therefore the problem can be reduced to the solution of one linear differential equation. 
We will then move to more interesting cases where, even in fixed numbers of dimensions, more than one 
master integral remain linearly independent and one cannot avoid the problem of solving higher order 
differential equations which give rise to more complicated mathematical structures. 


4 Explicit examples 

In the previous section we outlined the main ideas behind this paper. We argued that the IBPs might 
degenerate in the limit of fixed (even) integer numbers of dimensions d —>■ 2 n, such that some of the master 
integrals become effectively linearly dependent from each other. While this is a very well known fact, we 
argued that this degeneracy, if present, can be used in order to simplify the system of differential equations 
satisfied by the master integrals. In this section we will present many explicit examples of this simple 
idea. We will start by studying the two-loop sunrise graph with one massive and two massless propagators, 
Section l4Tl and a two-loop triangle with three off-shell legs. Section IT^ In both examples there are only 
two master integrals and by studying the IBPs in fixed even numbers of dimensions, one relation can be 
found, allowing to decouple the differential equations in that limit. We will then consider the case of the 
two-loop massive sunrise. Section H751 and of a non-planar two-loop triangle. Section HlTl In both cases not 
all equations can be decoupled, and a minimal bulk of two differential equations remains coupled, giving rise 
to elliptic functions. We will then study the case of a two-loop massive triangle with three master integrals. 
Section 14.51 Here, similarly to the non-planar two-loop triangle, there are three master integrals. In this 
case, nevertheless, the differential equations can be completely decoupled and the solution can be written 
in terms of MPLs. Finally, as a last example, we will move to the three-loop massive banana graph IT751 
In this case, we will study different possible mass-arrangements of increasing complexity, showing how the 
number of master integrals changes consequently, and how our method allows to determine easily which 
subset of master integrals can be immediately decoupled from the differential equations. 


4.1 The two-loop sunrise with one massive propagator 

Let us start off by considering the case of the two-loop sunrise with one massive and two massless propa¬ 
gators. We define the following set of integrals belonging to its Feynman graph 


/(d;ni,n2,n3,n4,n5) 





{k-pY^{i-pY^ 


m?) 


2'v"3 


where p^ = s \s the momentum transfer. The integration measure is defined as 


(4.1) 


= C{d) 
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(27r)'^ ’ 


(4.2) 





and the explicit form of the function C{d) is not relevant for the considerations below. Note that this 
Feynman graph does not contain any sub-topology. We keep explicit only the dependence on the space- 
time dimensions d and on the powers of the denominators and scalar products, which will be important 
for what follows. Performing a usual reduction through IBPs one finds two independent Mis, which can be 
chosen as 

Xi(d;s) =/(d; 1,1,1,0,0), l 2 (d; s) =/(d; 1,1, 2,0, 0). (4.3) 

Using the methods outlined in the previous sections we can now derive the differential equations fulfilled 
by Ii and I 2 in the momentum squared s. This step can be performed automatically using, for example, 
Reduze 2 [9], and we end up with the following 2x2 linear system 

dh 
ds 
dl2 
ds 

As one can immediately see, the equations are coupled for any even value of the dimensions c@. 

It is well known that these integrals can be computed as a series expansion in d —>■ 4 in terms of HPLs 
only, see for example |47j . Let us then try and use the ideas outlined in Section |3] in order to decouple the 
differential equations in the limit d —>■ 4. First of all note that, in the limit d —>■ 4, both master integrals 
are UV divergent and in particular they both develop a double pole 

+ (d^ (4; s) + 0(1) (4.5) 

^2 (d; S) = 4-"^ (4; s) + 4 -') (4; s) + 0(1). (4.6) 

Moreover it is easy to show that any integral of the form (SID can at most develop a double pole 
in (d — 4). Equipped with these consideration, let us now produce the IBPs for this Feynman graph as 
described above but, instead of solving them keeping the full dependence on the parameter d, we can set 
d = 40 As we discussed in detail in the previous section, this is equivalent to expanding the IBPs in 
Laurent series, and considering the first of the chained systems of equations obtained, namely the one 
corresponding to the double pole in (d — 4). Following the arguments of the previous section, we would 
expect to find a degeneracy of the two master integrals in d = 4, which should then allow us to decouple 
the two differential equations. As expected the two masters (|4.3I1 become linearly dependent 

4-^'>(4;s) = ^4-^H4;s). (4.7) 

As discussed in Section [31 such a relation must come from a corresponding d-dimensional IBP. Indeed, if 
one considers the original d-dimensional system of IBPs and solves it for the two masters, it is easy to find 
the following relation 

/(d;2,l,I,0,0) = -^^^^((3d-8)Ii(d;s)-4m2 2:2(d;s)) . (4.8) 

\d — 4J s — 

In the limit d —>■ 4 Eq. (14.81) trivially generates Eq. (14.71) . In the notation of Section [3] we can write this 
relation as 

(3d-8)Ii(d;s) - 4 m 2 l 2 (d;s) = C>(d-4), (4.9) 

or, equivalently, keeping also in the right-hand side only terms of 0{d — 4), 

Xi(d; s) — m^X 2 (d; s) = 0(d — 4), (4-10) 

®On tlie otlier hand, the equations become triangular as d —^ 3. 

^®The possibility of solving IBPs for fixed values of d is already implemented in the development version of Reduze 2. 


(d - 3) m2 

- M - -L2 

s s 

(d-3)(3d-8) /I 
2m2 Is 


s — 


Xi + 


2(d-3) (3d- 8 ) 


s — m^ 


2 s 


X2. 


(4.4) 


9 














recalling that this does not mean that this linear combination is of order 0{d — 4), but that it becomes 
zero if we fix d = 4 in the IBPs. 

In this particular case, since the Feynman graph under consideration does not have any sub-topologies, 
Eq. (14.71) can be seen as a real relation between the highest poles of the two master integrals. This relation, 
which is naturally derived from the IBPs only, can be easily verified by computing the highest poles of the 
two master integrals. A very simple exercise gives 




l2{d;s) 


I 

(d-4)2 

I 

(d-4)2 



((d-4) 

Jd^) 


(4.11) 


in agreement with Eq. dUZl). The overall normalisation of Eq. (I4.11|) is of course arbitrary and it has to 
do with the choice for the integration measure (14.21) . Let us now try and exploit this relation in order to 
simplify the system of differential equations (14.41) . We perform the change of basis from the “standard” 
Mis Ii{d; s), l 2 {d; s), to the new Mis defined as 


Ji{d]s) =Xi{d]s) -m^ I 2 (d; s), 7 ^ 2 (d; s) = Xi(d; s). 


(4.12) 


Note that in this case the first of the two masters in (14.121) has only a single pole in (d — 4) due to the 
exactness of relation (I4J1) . As second master integral we can choose any of the two and here we performed 
simply a random choice picking Ii (d; s). Choosing I 2 (d; s) would indeed lead to equivalent results. Deriving 
the differential equations for the new basis we find 


dJi 

ds 



i + (d-4) 
s 




dj2 

ds 


+ (d-4) 



3 


_2(s- 

(d-4) 


m?) 

■J2- 


S 


i + |(d-4) 


5 — 



J 2 


(4.13) 


Equations ()4.13|) confirm the discussion in Section [3] and can therefore be seen as of the main result of this 
paper. Let us have a closer look at these two equations and compare them to ()4.4|) . We note immediately 
that the equations are not in canonical form. On the other hand, the matrix of the system does become 
triangular in the limit d —>■ 4, and in particular the master ^2 appears in the differential equation for 
integral 77i multiplied by an explicit factor (d — 4), as predicted in p.8l) . For any practical purposes this 
is enough, since it means that one can expand the differential equations as Laurent series in (d — 4) and, 
order by order, first solve the differential equation for J7i by simple quadrature, and then use this result as 
input for the differential equation for J72, which can in turn be solved by quadrature. Needless to say, this 
procedure can in principle be iterated up to any order in (d — 4). 

A comment is in order. In this simple example the relation found by studying the IBPs in d = 4 can be 
interpreted as an actual relation between the double poles of the two master integrals. Very often the first 
poles of arbitrarily chosen Mis are either constants or very simple rational functions. One might therefore 
naively think that, by evaluating explicitly the poles of a given set of Mis, one could simply look for simple 
relations among the latter. Such relations would indeed contain only simple rational functions. It is well 
known, though, that in several cases the poles of a master integral can be represented entirely through 
its sub-topologies. If this is the case, a relation between the poles of the masters would be useless, as it 
would not bring any new information as far as the master integrals are concerned. In order to achieve 
a decoupling one must therefore use a relation which is contained in the IBPs and as such represents an 
effective degeneracy of the master integrals in the limit d —>■ 2 n, with n £ N. 


4.1.1 Simplification of the differential equations in d = 2 

It is interesting to see what happens by repeating the same analysis for the graph (sm in d = 2 instead 
of d = 4. Again, an easy analysis of the two master integrals ()4.3I) shows that they both develop a double 
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pole in d = 2, which in this case is of IR origin 


= Jd^ + (d^ (2; s) + 0(1) (4.14) 

= (d^ + (d^ X^-i) (2; s) + 0(1). (4.15) 

As before, one can easily see that also in this case all integrals of the form (|4.1D can develop at most a 

double pole in (d — 2). We can proceed and generate the IBPs for generic d and then expand them as 

Laurent series, this time in (d — 2), starting from l/(d — 2)^. We are then left with a series of chained 
systems of IBPs, each for a different order in (d — 2). As in the previous case, we can now focus on solving 
the hrst system, corresponding to the double pole. This again is equivalent to considering the original 
system of IBPs and simply fixing d = 2. Upon doing this one immediately sees that once more the two 
Mis become linearly dependent 

l(-'^(2;s) = ^-^x(-')(2;s). (4.16) 

s — 


As for the previous case, relation (14.161) must come from a corresponding d-dimensional relation. Indeed, 
if one solves the IBPs in d dimensions one finds, among the others, the following relation 


/(d; 1,1,1,1,0) 





It is clear that in the limit d —>■ 2 Eq. (14.171) generates instead Eq. (I4.16|) . Proceeding as above, we can 
choose as new basis of Mis 


Ji(d;s) = Xi(d;s) - (s - m^)X 2 (d;s), J 2 {d]s) = Xi(d;s). 

Deriving the differential equations for J\ and Ji one finds immediately 


(4.17) 


^ = (d-2) 

2 

3 

s — 2 s 

dj2 

1 

1' 

Ji + 

ds 

s — m'‘ 

s 


Ji + (d - 2) 
(d - 2) 1 


3(d- 2) 

2 ^ 

J 2 . 


J 2 


(4.18) 


Again, as expected from the arguments of Section [31 we see that the differential equation for Ji decouples 
from the one for S 2 in the limit d —>■ 2, respecting the same pattern described in equation (13.8|) . Once more 
for every practical purposes this is enough to reduce the solution of the system of differential equations to 
iterated integrations by quadrature. 


4.2 A two-loop triangle with three legs off-shell 

Let us consider now a massless two-loop three-point function with three legs off-shell. The problem has 
been widely studied in the literature, mainly in the context of vector boson pair production |48H51j , and it 
is well known that this Feynman graph can be reduced to two independent Mis, which can be integrated 
in terms of MPLs only. We define the Feynman graph as follows 


/(d;ni,n2,n3,n4,n5,n6,n7) = 


= / 



((fc _ ;)2)"2 _ p^)2)"4 , 


(4.19) 
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where pf = mf, p\ = rn^ and = {pi +^ 2 )^ = s. 

We used Reduze 2 in order to reduce this graph to two independent Mis 


= /(d; 1 , 1 , 1 , 1 , 0 , 0 , 0 ), l 2 (d; s, mljiri^) = /(d; 2 , 1 , 1 , 1 , 0 , 0 , 0 ). 


(4.20) 


We can then proceed and derive the differential equations for these two Mis. As always we neglect the 
sub-topologies throughout. In this particular case the latter are simple two-loop corrections to massless 
two-point functions which have been known analytically for a very long time. 

The homogeneous part of the differential equations in the momentum transfer s reads 




2 2 ^ <^^1 (d — 4)(mf — -I-((3d — 8 ) 771 ^ — 3(d — 4 )to|) s-I-2(d — 4)s^ 2 


2 s 


Ii -f 2 s X 2 

(4.21) 


2 2 ^dl 2 (10-3d)((d-3)(m?-TOi)-k(2d-7)s) 

P(s,mi,m2) -j— = -^--- -Ii 

as Zs 

(d — 6)(mf — m^)^ + ((22 — 5d)mi -I- (d -f 2 ) 1112 ) ^ ~ 2(d — 2) s^ 

2^ 


X 2 


(4.22) 


where we defined the polynomial 


P{s, ml, m 2 ) = m^P {s — ml)^ — 2m^(s -I- ). 


(4.23) 


The equations are coupled in the limit d —>■ 4. Again, as for the sunrise studied in Section 14.11 the 
integrals can develop at most a double pole in (d —4). Instead of performing a complete Laurent expansion 
of the IBPs, we generate them and then fix explicitly d = 4 before solving them. This is enough to check 
whether the two Mis degenerate in this limit. By solving the IBPs one finds that this is precisely the case 
and the following relation is extracted 

Ii(d; s,ml,ml) + sl 2 (d; s, m\,m\) = C>(d — 4), (4.24) 

where again we used the notation introduced above, indicating that the combination (14.241) becomes zero 
if we fix d = 4 in the IBPs. Of course, also in this case, if we had expanded the IBPs as Laurent series 
starting from the double pole, relation (14.241) could have been interpreted as a relation between the double 
poles of the two master integrals 

T[~'^\x,s,m\,m\) + sl2~'^\d; s,ml,ml) = 0. (4-25) 

Note, nevertheless, that this time the relation is not exact, differently from (IQ) . since the sub-topologies 
might in general contribute modifying (I4.24|) . Eq. (14.241) is anyway sufficient to decouple the homogeneous 
part of the differential equations. We proceed as above and define the new basis 

y7i(d; s, mf, m 2 ) = Ii(d; s, m\,m\) + sl 2 (d; s, m\,ml), J 2 {d', s, m^, m^) = Ii(d; s, m\,ml). (4.26) 


Deriving the differential equations satisfied by (j4.26l) we find 


P{ 


s, m,, m 


N dJi 


Pi 


s, mr, m: 


> dJ2 


(d — 4)(mi — ml)^ + ((22 — bd)m\ -|- (d — 2)m|) s — 2(d — 3)s^ 
-2^- 

(d — 4) (3(d — 5)ml + (11 — 3d)ml — 3(7 — 2d)s) 
-2- 

o™2 ^ , {d-4:){s + ml-ml){2s + ml-ml) ^ 

^ m^ Ji --- J2 ■ 


(4.27) 


Again, as expected, we see that the differential equation for Ji contains an explicit factor (d — 4) 
multiplying the second integral J 2 - The result is consistent with the general structure described in Section[3l 
Once more we want to stress that, as expected, in this case all Mis can be integrated in terms of MPLs 
only. 
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4.3 The two-loop massive sunrise 

In the previous sections we considered two simple examples of 2 x 2 systems where the two equations 
could be decoupled in the limit d ^ 4, such that the problem could always be reduced to integrations by 
quadrature. As we already discussed this is not always possible and the first known case where at least two 
differential equations remain coupled is the two-loop massive sunrise. The two-loop massive sunrise graph 
is defined as follows 


712,713,714, ns) 



{k-pY^{l-pY^ 

— (P — 7n|)”^ ((fc — I +p)‘^ 



n3 ■ 


(4.28) 


This integral has been studied widely in the literature and in particular a lot of attention has been devoted 
to the differential equations that it fulfils. In the general case where all three masses assume different 
values, a normal reduction through IBPs shows that all integrals can be expressed as linear combinations 
of 4 independent Mis, which can be chosen to be 


Ii(d;s) =/(fi;l,l,I,0,0), 
Isid; s) = /(d; I, 2,1, 0,0), 


l2{d\s) =/(d;2,I,l,0, 0), 
l 4 (d; s) = I{d; 1,1, 2,0, 0). 


(4.29) 


In [TS] it was shown that these integrals fulfil a coupled system of 4 linear first order differential equations 
in d dimensions. The system remains coupled in the limits d —>■ 2 ti, where ti € N. It was lately shown 
in [52], using algebraic geometry methods (and as such a priori orthogonal to the IBPs), that the scalar 
integral Ji(d;s) satisfies a second-order Picard-Fuchs differential equation in d = 2. This suggested the 
possibility of finding a proper change of basis of Mis, in the sense of the IBPs, such that two of the four 
differential equations satisfied by the latter would decouple in the limit d ^ 2. Since the four Mis in (14.291) 
are finite in d = 2, it appeared natural to try and obtain the decoupling of the differential equations by 
finding new relations among the Mis, valid strictly only for d = 2. Such relations can be found using the 
so-called Schouten Identities introduced in m- In that reference the Schouten identities are introduced 
and the case of the two-loop sunrise with different masses is worked out in detail. It is shown that, as 
expected, in d = 2 only two master integrals are linearly independent. This allowed to recover the second 
order differential equation found in reference |52] in a completely independent manner. In this section we 
will show that those results can be even more easily re-obtained using the methods described in this paper, 
and namely by solving the IBPs for the massive sunrise in d = 2. The Schouten identities can be imagined 
as a tool for extracting this piece of information from the IBPs and are, in this respect, equivalent to the 
study of the IBPs in fixed number of dimensions. We will show another example of this equivalence in 
Appendix 

Since the algebra in this case is rather heavy due to the large number of scales, we will only report the 
result of the solution of the IBPs in d = 2, referring to [41] for their use to simplify the system of differential 
equations. As already discussed above, solving the system with d = 2 is in general easier and, as expected, 
two of the four Mis degenerate, leaving only two linearly independent Mis. By choosing as Mis Ti(2; s) 
and 22 ( 2 ; s), we find the following additional relations (as everywhere else we neglect the sub-topologies for 
simplicity) 


m\P{s,m\,m\, m\)I^{2-, s) = + m\ — m\ — s)Xi(2; s) 

-I- 777^ (^mf — 37772 -|- 2777^(7772 — 777^ — s) -|- (m§ — s)^ -|- 27772(7773 -|- s)) 22(2; s) 

7773 P(s, 777J, 7772, 7773)24(2; s) = (777^ — 7773)(777^ — 777^ -|- 7773 — s)2l(2; s) 

-I- nil (™i + ~ 37773 -I- 27772(7773 — s) -I- 27773 S -I- — 2m\{m\ — m\ + s)) 22(2; s ), (4.30) 
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where we defined the polynomial 


P(s, m\,m 2 , TO 3 ) = (—3mf + + (TO 3 — s)^ — 2m2(m^ + s) + 2m\{m\ + TO 3 + s)). 

As we discussed in Section [3J we expect such relations to come by d-dimensional IBPs with an overall 
factor l/(d — 2). Indeed by studying the reduction to Mis in d dimensions it is easy to find the following 
two relations 


0{d — 2) = [(d — 3){2mi — ml — m^) — (d — 2)s] Ii(d; s) -I- 2 ml(s — m\)l 2 {d\ s) 

+ ml{—3m\ +ml + 3m\ — s)l 3 (d; s) -I- + 3ml + m\ — s)Ii{d; s)| 

0{d — 2) = -| [(d — 3){ml — 2mf -I- m 3 ) — (d — 2)s] Ii(d; s) — 2 ml{s — ml)Iz{d] s) 

-I- m\{—m\ + 3ml — 3m\ + s)l 2 (d; s) + + 3m\ — m\ -1- s)X 4 (d; s)|. (4-31) 

Relations (14.311) can be compared with the corresponding formulas (3.14) and (3.15) of [H]. It is easy to 

see that they are identical in the limit d ^ 2 , the only difference being the absence of the terms coming 
from the sub-topologies, which we are neglecting here. These two relations (and in particular their limiting 
value as d —2) can be used, as described in Section [3l in order to decouple two of the four differential 
equations of the two-loop massive sunrise graph, by choosing as new basis of master integrals 

Ji(d; s) = Ii(d; s), J 2 {d] s) = X 2 {d] s) 

Js{d-, s) = —{2ml “ ”^2 “ ™ 3 ) 2 ii(d; s) + 2ml{s — ml)l 2 {d; s) 

+ ml{—3ml -I- m^ -I- dmg — s)Xz{d] s) + ml(—3ml + 3ml + mg — s)l 4 (d; s) 

f 74 (d; s) = —(mf — 2mf -I- m 3 )Ii(d; s) — 2mf(s — m 2 )X 3 (d; s) 

-I- ml(—ml + 3ml — 3ml + s)X 2 {d; s) + ml{—3ml + 3m\ — mg -|- s)l 4 (d; s). (4.32) 

We do not give the explicit form of the differential equations, referring to [41] for further details. In 
comparing, note that the basis presented here differs from the one in [JT] by the absence of sub-topologies 
and by orders 0{d—2). Furthermore we want to stress, in relation to the discussion in Section|3J that using 
this basis produces an overall factor (d— 2 ) in front of the two differential equations for J^ld] x) and J^ld] x). 
This implies that one has, at every order in (d — 2), only one second order differential equation (needed 
to solve the block of Ji{d]x) and J 2 {d\x)), plus two integrations by quadrature (required to determine 
J‘i{d]x) and Ji{d-,x)). 

4.4 A two-loop non-planar crossed vertex 

As a further application, let us consider a two-loop non-planar crossed vertex with two massive propagators. 
This graph is topologically completely unrelated to the two-loop sunrise and was studied thoroughly in . 
There it was shown that it can be reduced to three Mis, which would therefore be expected to satisfy a 
system of three coupled differential equations. In [53] a basis of Mis was found such that one of the three 
differential equations decouples from the other two in the limit d —>■ 4. This reduced effectively the problem 
to that of solving, for every order in (d — 4), a second order differential equation, plus an integration by 
quadrature for the third MI. In this section we would like to study this Feynman graph with our method 
and show that the decoupling found in |53| comes as well from a degeneracy of the master integrals in d = 4 
which can be read off directly from the IBPs. Following [53] we define the Feynman graph as follows 
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P 2 


= / D‘^kD'^1 


[k-p2Y 


(P _ m^)ni {12 _ {{k _ p^)2)"3 {{I _ {{k-l- ’ 

(4.33) 


with pj = 4*2 = 0 (pi + p 2 )^ = s. It is easy to verify that this topology can be reduced to 3 Mis, for 
example 


Ii(d;s) =/(d; 1,1,1,1,1,1,0), l 2 (d;s) =/(d;2,l,l,l,l,l,0), JsCd; s) =/(d; 1,1, 2,1,1,1,0). (4.34) 


Let us derive the differential equations in the momentum transfer s, neglecting as everywhere else all 
sub-topologies. We get 


dll {d-&) 2m^ 

— — - -Ll - -L2 , 

as s s 

dl 2 {d — 5){2d — 9)(s — 4m^) ^ 14(d — 4)m"^ — (5d — 13)m^s — 2s^ ^ 

ds s{s — m‘^){s + 8m'^) ^ s{s — m'^){s + Sm'^) ^ 

2 (d - 4) 
s(s-|-8m^) 

dis 
d s 


Looking at these equations we see immediately that I 3 is already decoupled from the other two in the 
limit d —>■ 4. This means that, with this basis, the problem is reduced to that of solving, at every order 
in (d — 4), a coupled system for Ii and 12- With the explicit solution for the latter, one can then obtain 
I 3 integrating its differential equation by quadrature. It would be then interesting to know whether this 
decoupling is also due to a degeneracy of the Mis in d = 4. Moreover it would be even more interesting to 
verify whether a new basis could be found, for which the differential equations become completely triangular 
as d —>■ 4, reducing even further the complexity of the problem. 

Following the recipe described above, we can try and solve the IBPs for this Feynman graph for d = 4. 
A word of caution is required here. The three Mis selected above have Laurent expansions in (d — 4) which 
start at different orders, in particular one can easily find (for example using sector decomposition [^) that 
the first two masters are finite, while the third develops a cubic pole 

Ii(d^4;s) = 0(l), l2{d^4;s) = 0{l), I^d ^ 4; s) = O (4.36) 

Nevertheless this poses no practical obstacle to the applicability of the method presented in this paper. 
As we already discussed in general, by expanding the system of IBPs in Laurent series around d = 4, we 
will get, at every order in (d — 4), an independent system of differential equations whose homogeneous 
part (i.e. the one containing the order of the Mis under consideration) has always the same form, while 
the non-homogeneous part will of course change and, in particular, depend on the previous orders of the 
expansion (and on the sub-topologies, which we neglect). What we are interested in is, indeed, only the 


2(d - 5)(2d - 9)(s -I- 2m^) 2m^ ((24 — 5d)m^ -|- (21 — 4d)s) 

^1 ”1- 7 -— TT , —9\-^2 


s(s — ■w?){s -I- Sm^) 
2(3d — 8)rrY + {d — 3)s 
s(s -I- 8 m^) 


s(s — ■rn?){s 8m‘^) 


X 3 . 


(4.35) 
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homogeneous part of this system. Fixing d = 4 is therefore enough in order to determine whether, for any 
order of the expansion, the Mis become linearly dependent. Upon doing this we find only one relation 
between the three masters which reads 

Ti{d]s) + {bm^ + s)X2{d]s)+Zm^l3{d]s) = 0{d - 4) , (4.37) 

where, as always, we mean that this combination becomes zero when we set d = 4 in the IBPs. 

Equivalently, one can also proceed in a more formal way, expanding all IBPs in Laurent series starting 
from the triple pole up to the finite piece, and supplementing the piece of information on the highest poles 
of the Mis (|4.36l) . Upon doing this, one obtains four chained systems of IBPs (one for every oder in (d —4)), 
which can be solved bottom-up starting from the one corresponding to the highest pole. Since the first 
two masters are finite, the first three systems give no information on the latter, while the fourth system 
(corresponding to the finite piece of the Mis) produces the relation 

lf\4- s) -h (5m2 -h s)J^°^(4; s) s) + 1,1,1,1,1,1,1) = 0 . (4.38) 

Indeed relations (14.371) and (14.381) are identical up to the presence of the previous order in the expansion 
of the integral /(d; 1,1,1,1,1,1,1). If we had solved the IBPs in d dimensions, this integral would have 
been of course expressed in terms of the three masters (14.351) . Solving the system in d = 4 instead does not 
allow to express this integral in terms of the other three, but this comes with no surprise and can be very 
well understood in terms of the degeneracy of the system of IBPs in this limit. 

By studying explicitly the integral /(d; 1,1,1,1,1,1,1) it is easy to see that it is also finite in d = 4, 
namely 

/(d^ 4; 1,1,1,1,1,1,1) = 0(1), ^ /(-i)(4;l,l,l,l,l,l,l) = 0. 

With this piece of information one recovers again relation (14.371) . which was found by simply solving the 
system of IBPs in d = 4. Since only one relation has been found, which moreover involves all three masters 
Ii, X 2 , and X 3 , we have no way to decouple more than one Mis from the system. What we mean here 
is that, since in system (j4.35l) only the differential equations for Xi and I 2 are coupled, if we had found 
one relation but involving only Xi and X 2 , we could have used it to decouple this block of the system. An 
example of this is given in Section 14.51 This is not the case and we therefore expect the minimal number 
of coupled integrals in d = 4 to be two, giving rise to a second order differential equation, as for the case 
of the two-loop massive sunrise, see Section 

As an exercise, we can try to change basis also in this case exploiting the piece of information found 
in (14.371) . We expect to end up with a new system of differential equations, where nevertheless again two 
out of three equations are coupled as d —>■ 4 (and as such practically equivalent to (14.351) 1. showing that 
the system cannot be further simplified. Let us introduce the new basis 

J7i(d; s) = Xi(d; s), J72(d; s) = X 2 (d; s), J 3 (d; s) = Xi(d; s)-|-(5 -I-s)X 2 (d; s)-I-3 X 3 (d; s). (4.39) 

Deriving the differential equations and neglecting all sub-topologies we get 

dJi (d - 6 ) 2 

- — - Jl - J 2 

as s s 

dJ 2 2(d — 4)(s — TO^)-I-3(d — 5)(2d — 9)(s — 4m^) 
ds 3s(s — m^)(s-I-8 to2) ^ 

52(d — 4)m^ — (23d — 7I)m^ s — 2(d — I)s^ 2 (d — 4) m? 

3s(s — m^)(s-I-8m2) J 2 + 


dJs (d-4)(6d-29) (d-4)(s-I0m^) (d - I) J 3 

ds Qm'^s ^ 'dm?s ^ 3 s 


(4.40) 
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which is again a system of three differential equation, two of which remain coupled in the limit d ^ 4, giving 
rise to a second order differential equation for one of the two coupled masters^ We note that the new 
system (|4.40p , compared with the previous one (14.351) , has a slightly different structure. As for the previous 
cases that we analysed, once we switch to the new basis defined through the IBPs degeneracy (14.371) . the 
differential equation for the new master develops an explicit factor (d — 4) in front of the other two 
masters Ji and J 72 , as predicted in equation (13.81) . 

4.5 A two-loop massive triangle with three master integrals 

Before moving to a three-loop example, let us try and see what happens in a case similar to the one 
studied above, i.e. a Feynman graph reduced to three master integrals, but where the system of differential 
equations can be completely triangularised as d 4. Let us consider the following two-loop massive 
triangle 


/(d; ni,n 2 , na, n4, ns, ne, n^) 



= / 


(fc-pir {i-pir (i-qY 


[P - ((fc - {{k - piY - {{k -Pi- P2Y - m?) 


2\n4 5 


(4.41) 


with two legs off-shell, namely = (pi + P2) = s, p\ = 0, pi = q^. This graph has been studied in the 
context of the QCD corrections to id —>■ Zj in [^[^. Similarly to our previous example it is reduced 
to three master integrals, such that we are dealing with a system of three differential equations. We start 
from an arbitrarily chosen basis of master integrals 


Xi(d;s,(?^) = /(d; 1 , 1 , 1 , 1 , 0 , 0 , 0 ), J 2 (d;s, g^) =/(d; 1 , 1 , 2 , 1 , 0 , 0 , 0 ) 


X3(rf;s,g2) = /(d; 1,1,1,2,0,0,0). 


(4.42) 


The masters depend on three variables s, and m^, and therefore on two independent ratios. For simplicity 
we will consider only the differential equations in s, while all considerations done here work identically for 
the differential equations in the other variables. In order to simplify as much as possible the formulas we 
write explicitly only the order zero of the homogeneous differential equations in (d — 4), which is also the 
bulk which we need to simplify. The equations read 


^ n- 1 ^ \ 1 


s — 5) 


d TTlP' 

as s[q^ ~ ~ 


X2(d; s) + 


— T (P- - g) 

OS s (s(g^ ~ -5) — m^q^) 


X2(d; s) -|- 


q^s — mr{ 2 q^ + s) 

(g2 — s)(s(g^ — s) — rrPq^) 
wP{s^ + sq^ — g^) 


s(g^ — s)(s(g^ — s) — m?q^) 


^ 3 ( 4 ; s) -)- 0 (d — 4) 
X 3 (d; s) + C)(d-4) 


(4.43) 


where the dependence from g^ is left as implicit in the master integrals for ease of notation. One can 
immediately see that only two of the differential equations are coupled. One should in principle first solve 
the 2 x 2 coupled system for X 2 (d; s) and X 3 (d; s), and then, with the latter as an input, one could attempt 
to solve the differential equation for Xi(d; s) by quadrature. 

Let us try now and study the IBPs in the limit d —> 4. By solving them as discussed in the previous 
sections one sees that the master integrals X 2 (d; s), X 3 (d; s) become linearly dependent in this limit and one 
finds the relation 


(g^ — s)(s — 2 m^)X 2 (d; s) -I- [s(s -|- 2m^) — 2 g^ (s — m^)] s) = 0 {d — 4). (4.44) 

1531 it was shown that the homogeneous part of the second order differential equation satisfied by the scalar master 
integral Xi(d; s) is equivalent to that of the two-loop massive sunrise with equal masses. 


17 




















As for the case of the non-planar triangle studied in Section 14.41 we find only one relation, while we 
have three master integrals. One of the three masters nevertheless is already decoupled, and moreover 
relation (|4.44l) involves only X 2 (d; s) and Xsld; s), which are precisely the two coupled integrals. We expect 
this therefore to be enough to decouple the system. We define the new basis 


=Ii(d;s), J 2 {d] s) = l 2 {d\ s) 

( 9 ^ - s)(« - 2 ^ 2 )^ , s{s + 2m'^)-2q^{s-m^)^ 

jTslct, 5) — . X2(a, s) -|- . X3(a, 5) . 


(4.45) 


where the 1/m^ has been added for dimensional reasons. Deriving the differential equations for this new 
basis, and keeping again only the first order in [d — 4) we find 


d ^ ^ . 

= 


— s 


Ji{s]d) 


s{q^ — 4771^) 


2 q^{s — m?) — s{s + ’m?) 


J 2 {d-, s) 


(2 m?{q^ + s) — sq^)\ 


s{q^ — s){2q^{s — m?) — s(s + im?)) 


Jsid', s) + 0{d — 4) 


N _ q'^{2m'^ + {s-2m‘^)s) _ 

ds ^ (q^(s — m^) — s^)(2 q^(s — m^) — s(s + 2m^)) ^ 

(s nn? — q^{s — rn?)) rn^ 

{q^ — s){q^{s — m?) — s^){2 q^{s — m?) — s(s + 2m^)) 


J 3 (d; s) + 0(d-4) 


^J3(d;s) = ^^-^J3(d;s) + 0(d-4). (4.46) 

os s[q^ — s) 

As expected the system of differential equations becomes triangular and, in particular, the equation 
for the new integral J 3 {d]s), defined through relation (14.441) . decouples from J 2 {d]s), following the usual 
pattern of equation (13.81) . One can then proceed, order by order in (d — 4), integrating by quadrature first 
the differential equation for J 73 (d;s), then the one for J 2 {d]s) and finally the one for Ji{d\s). As a last 
comment we want to stress that, if we had considered the system in d/dq'^, the same change of basis (14.451) 
would have indeed been sufficient to triangularise this one as well. 


4.6 The three-loop massive banana graph 

As last example let us consider a more complicated three-loop graph. We choose the three-loop massive 
banana graph, which is the natural three-loop generalisation of the two-loop massive sunrise. In the most 
general case this Feynman graph depends on the momentum squared = s and on four different masses 
mi, m 2 , m 3 and mi 

hid; ni, 712, n3,n4, ng, ng, 777 , ns, ng) = 

7)<^h _ (fci ■p)^^{k2 ■p)"^(fc3 ■p)"^(fci ■ fc 2 )""(fci ■ fcs)"" _ 

^ ^ ^ (fc^ — mi)"i(fc| — mf)"2(fc| — 77i2)"-3((fcj^-(_/c2 + A:3 — p)2 — mi)"^ ’ 

(4.47) 

where the subscript 4 indicates that the four masses are all different. In the two-loop case there are 4 Mis 
when the 3 masses have all different values, which in turn degenerate to 2 Mis in the case of equal masses. 
On the other hand we saw that, irrespective of the explicit values of the internal masses, one is always left 
with only two independent Mis in d = 2 (14.301) . This allowed us to decouple two of the four Mis from 




18 



















the differential equations in the limit d —> 2 and prove that the scalar amplitude satisfies a second order 
differential equation in this limit. 

It would therefore be interesting to verify whether a similar behaviour can also be seen in the three-loop 
banana graph. Since in the general case with four different masses the algebra becomes very cumbersome, 
we will consider different cases of increasing complexity, namely increasing at every step the number of 
different internal masses and check how many Mis are found in d dimensions and how many can be decoupled 
in the limit d —>■ 2 . 


4.6.1 The equal-mass case 


Let us start considering the equal-mass case. We use the following notation 


/i(d; ni, 712,713, «4,715, «6, nr, ns, ng) 


hid; Til, 772,773, 774, 775, ^^6, 777, 778, 779) 


7714 = 7713 =7712=1711 =771 


where the subscript “1” indicates now that all masses have the same value. Running a reduction to Mis 
with a code of choice it is easy to check that there are 3 independent Mis in d dimensions which can be 
chosen to be 


Xi(d; s) = hid; 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 ), l 2 (d; s) = h (d; 2 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 ), 

X 3 (d;s) = Ji(d;3,I,l,I,0,0,0,0,0). 


The differential equations in the momentum transfer for these three Mis read 
dXi 3d — 8 ^ 4 777^ ^ 

=-Xi-X 2 

ds 2 s s 


dX 2 (3d-8)(2d-5) (d - 4 ) s - 8 ( 2 d - 5)7772 1 

~J Q 2 C 2 O 

ds 8 s 777^ 8 S 777 ^ 2 


(4.48) 


dXs _ (2d - 5)(3d - 8) (I6(lld - 37)777^ -f 4(32 - 9d) s -h (d - 4) s^) 
ds 32 777^ s (s — 47772 )(s — I 67772 ) ^ 

[64 (440 -h (47d - 289)d) - 16 (668 -h d(62d - 409)) m^s -f 16(d - 4)(4d - 13 ) 7772 s^ 

32 m? s (s — 47772 )(s — I 67772 ) 

[I024(d - 4 ) 777 ® -H 192(27 - 8d)777®s -h 96(2d - 7)7774s2 - 4(d - 4 ) 7772 s®] 

32 7772 S (s — 47772 )(s — I 67772 ) 


(d-4)V]^ 
- 2-2 

(4.49) 


and we can easily verify that, in spite of the fact that X 3 does not appear in the first equation, the system 
is still coupled as d 2. Trying to solve the system of IBPs in d = 2 shows no further degeneracy and 
therefore we conclude that the system cannot be further simplified with our method. Having a system 
of three coupled first-order equation means that we can rephrase it as a third-order differential equation 
for any of the three masters, and in particular for the scalar amplitude Xi(d;s). The fact that the scalar 
amplitude fulfils a third-order differential equation is in agreement with the findings in [52) . Deriving the 
third-order differential equation satisfied by Xi(d; s) we find 


Dfh{d;s)=Q, (4.50) 

( 3 ) 

where the d-dimensional third order differential operator ' reads 
( 3 ) d® 3 (64777"^ -I- 10(d — b)rn?s — (d — 4)s2) d2 

ds® s(s — 47772 )(s — 16777^) d s'^ 

(d — 4)(lld — 36) s 2 — 64(d — 4)d 777'^ — 4 (216 -I- d(7d — 88)) m? s d 
4 s2(s — 47772 )(s — 167772 ) ds 

(3-d)(3d-8) ( 2 (d-h 2 )?T 72 -h(d- 4 )s) 

4 s 2 (s-47772)(S-I67772) ’ ^ 
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and all sub-topologies are neglected as always. In the limit d —>■ 2 the differential operator simplifies to 


(3) d? 6 (s^ — ISm^s -I- 327^"^) d? (7s^ — 68m^s -|- 64 to'^) d 1 

^ ds^ s(s — 4m^)(s — 16m^) ds^ s^(s — 4m^)(s — 16m^) ds (s — 16m^) ’ 

which is in agreement with |52) . 


4.6.2 The case of two different masses 

Let us move now to a slightly more general case and let the masses take two different values. There are 
two possible arrangements, which we call and defined as follows 


I^(d; ni, 712, ns, 774,775,775,777, ns, ng) = hid; 77i, 772, ns, 774,775, ng, nr, ns, ng) 
I2 [d; 77i, 772, ns, 774,775,776, nr, ns, ng) = I/^{d; tii, 772,773,774,775,775, nr, ns, ng) 


m 3 —m 2 —mi—ma ,7714—rrn, 


m2—7Tli—ma ,7714 — 7713—7711, 


The two configurations are intrinsically different and it makes sense to look at the two cases separately. 


A) In configuration A a reduction to Mis for generic d gives 5 independent Mis which can be chosen as 
If id; s) = If id; 1,1,1,1,0, 0,0, 0,0), if id; s) = if id; 2,1,1,1, 0,0, 0,0, 0), 


If id; s) = If id; 1,1,1,2,0, 0,0, 0,0), if id; s) = if id; 3,1,1,1, 0,0, 0,0, 0) 
lfid;s) = lfid;2,2,l,l,0,0,0,0,0). 


(4.53) 


B) In configuration B we find instead 6 independent Mis for generic d 


If(d;s) = /f(d; 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 ), 
If id; s) = If id;l,1,2,1,0,0,0,0,0), 
If id; s) = If id;2,2,1,1,0,0,0,0,0), 


If id; s) = If id;2,1,1,1,0,0,0,0,0), 
If id; s) = If id; 3,1,1,1, 0,0, 0,0, 0), 
l6^(d;s)=/f(d;2,l,2,l,0,0,0,0,0). 


(4.54) 


A natural question at this point would be how many Mis degenerate in the two cases in the limit d —>■ 2, 
and therefore what is the order of the differential equation satisfied by the scalar amplitudes if id; s) and 
if id; s) respectively. A naive expectation, based on the two-loop sunrise, would be to see in cases A and 
B, 2 and 3 Mis decouple respectively, such that the problem would reduce to the solution of a third-order 
differential equation, as in the equal-mass case. Unfortunately this naive expectation is not satisfied and 
we find that, by solving the IBPs in d = 2, in both cases four Mis remain independent, corresponding in 
principle to a fourth-order differential equation for the scalar amplitude in both mass-configurations. On 
the other hand, it is interesting to see that in both configurations, in spite of the different number of Mis 
in d dimensions, the problem can be reduced to an equation of the same order (i.e. four) in d = 2. 

Neglecting the sub-topologies we find in configuration A the following relation which allows to express 
the fifth master integral in terms of the previous four 


mlis-5ml-\-nil)lfi2;s) 


3 ml + nil-s , hlmi -b (777^ - s)^ - 67772(7772 -b 2s) 

mliml - s ) , 21777^ -b (777^ - s)2 - 67772(7772 -b s) 

---5-X3 Sj -I --X4 [Z, S) . 

6777 J 6 


■s) 

(4.55) 


In configuration B, instead, there are two different relations, which can be used to two express if and 
if in terms of the other four Mis in d = 2. We do not report the explicit solution of the IBPs in d = 2 
which looks rather cumbersome. As for the case of the two-loop sunrise, these identities originate from 
d-dimensional IBPs which degenerate in the limit d —>■ 2 due to an overall factor l/(d — 2). There are many 
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of these relations, but only two of them are linearly independent in the limit d —>■ 2, and they read (keeping 
only the first order in (d — 2)) 

0{d - 2) = [2(m^ + ml) - s] Ii(d; s) + [4to^(5to^ + 4to^) - 4 (2m^ + ml)s + s^] l 2 (d; s) 

+ Aml{2ml + ml - s)l 3 {d-, s) + 2771^ [4(m^ + ml){2ml - s) + s^] l4(d; s) 

+ Am‘l[2{ml+ml)-s] l5{d; s) + 8mlml{dml - s) Ie{d; s), (4.56) 


0{d — 2) = (—2m^ + 6ml — 3s) Ji(d; s) + [—20m^ + 8m^(7m^ — 2 s) + 3 s(s — 4m^)] l 2 (d; s) 

+ 12ml{ml — s)l 3 (d; s) + 2 m^ [— 8 m^(m^ — iml) — 4(m^ + 3m^) s + 3s^] Ii{d\ s) 

— 4m^(2m^ — + s)l 5 {d; s) + 32 ml ml{ml — s)l 6 (d; s). (4-57) 

We stress again that relations (I4.55p . (I4.56|) and (14.571) are not exact since all sub-topologies have been 
neglected throughout. These relations can be nevertheless used in order to derive new systems of differential 
equations where, for both A and B configurations, only 4 equations remain coupled in the limit d —^ 2. For 
example, in the case of configuration A we can take as new basis 

tA,^. ^Ar^. jA^^. ^ j-A^^. ^ 


J^^{d- s) = It{d- s) , j/(d; s) = I^(d; s), J^{d- s) = I^{d; s ), 


plus the new master defined as 


J 5 ^id] s) = ml{s - hml + mg)I^(2; s) - + "^6 -f x^{2-, s) 

12 mz 


51m4 + {ml - s )2 - &ml{ml + 2 s) ^ 


12 


1 ^( 2 ; s) 


mliml - s) ^ 21 m 4 -h (mg _ s)2 _ 6 to2(^2 

- a - 9 - ^3 1 ^) ■SJ ^- -^4 1^1 SJ . 

D m* 0 

Upon doing this one finds that the differential equation for the new master assumes the form 
= (d - 2) [c5i(d;s)Ji + C52(d; s)J2 + C53(d;s)Ji + C54(d; s)J2] + C55{d;s)J5 , 


(4.58) 


(4.59) 


where the functions Cij{d',s) are rational functions for the dimension d, the momentum s and the two 
masses, and are finite as d —)• 2. This insures, thanks to the overall coefficients (d— 2), that the differential 
equation for J 5 decouples completely from the other four, as expected. 

As far as configuration B is concerned, in order to achieve the complete decoupling of two out of the 
six equations, we can take as basis 

A _ ^^ 7 . 7 . A jf {d; s) = {d; s) , (4.60) 


Jf (d; s) = If (d; s), (d; s) = if (d; s), jf (d; s) = (d; s), 


together with 

Ji^{d;s) = [2{ml + ml) - s] Ii(d;s) + [Aml{5ml + dml) - A {2ml + ml)s + s^] 2 l 2 (d;s) 

+ 4m^(2m^ + ml - s)l 3 (d; s) + 2ml [4(m^ + ml){2ml - s) + s^] l 4 (d; s) 

+ 4m^ [2{ml + ml)-s] l 5 {d-, s) + 8mlml{Aml - s)lQ{d; s), (4.61) 


jf (d; s) = {-2ml + - 3s)Ii(d; s) + [-2Qml -I- 8ml{7ml - 2 s) -|- 3s(s - 4mb)] l 2 (d; s) 

+ 12 ml{ml — s)l 3 (d; s) + 2ml \—8ml{ml — 3ml) ~ 4(»Ba + 3ml) s + 3s^] l 4 (d; s) 

— 4m^(2TO^ — 6ml + s)l 5 (d; s) -I- 32mlml{ml — s)Ie{d] s). (4.62) 

Using this basis one obtains a new system of differential equations, where the two equations for and 
J7g^ develop an explicit overall factor (d — 2), such that, at every oder in the Laurent expansion, they can 
be solved trivially by quadrature. Order by order, once the result for the latter is known, one is left with a 
system of four coupled differential equations for the remaining master integrals. For compactness we prefer 
not to give here explicitly the systems of differential equations. 
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4.6.3 The case of three different masses 

Generalising even further we can check what happens if three out of the four masses are allowed to take 
different values. In this case there is of course only one possibility, which we choose to be 


/ 3 (d; ni, 712,713,714,715, Tie, 717 , ns, Tig) = h{d-, tii, 712 , 713 , 714 , 715 , rie, 717 , ns, ng) 


We start, as always, performing a reduction for generic d. The complexity increases and we find 8 
independent Mis 


Ii(d;s) =/ 3 (d; 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 ), 
l3(d;s)=/3(rf; 1,2,1,1,0,0,0,0,0), 
X 5 (d;s) =/ 3 (d; 3,1,1,1,0,0,0,0,0), 
l7(d;s) =/3(d; 2,1,2,1,0,0,0,0,0), 


l 2 (d;s) =/ 3 (d; 2 ,l,l,l, 0 , 0 , 0 , 0 , 0 ), 
X 4 (d;s) =/ 3 (d;l,l, 2 ,l, 0 , 0 , 0 , 0 , 0 ), 
X 6 (d;s) =/ 3 (d; 2 , 2 ,l,l, 0 , 0 , 0 , 0 , 0 ), 
X 8 (d;s) = J 3 (d;l, 2 , 2 ,l, 0 , 0 , 0 , 0 , 0 ). 


(4.63) 


We can then consider the system of IBPs for d = 2. It is easy to check that in this case 3 Mis 
degenerate, and therefore only 5 Mis remain linearly independent. We do not report here the equivalent 
of relations (14.551) (14.561) and (14.571) , since they are considerable more lengthy, but one can easily work out 
the reduction in d = 2 and find that, for example, X 6 ( 2 ;s), X 7 ( 2 ;s) and X 8 ( 2 ;s) can be written as linear 
combinations of the Xi(2; s),...,X 5 ( 2 ; s). Using the methods described above, 3 out of the 8 differential 
equations for this particular mass configuration can be decoupled in the limit d —> 2 , and one can in 
principle derive a fifth-order differential equation for any of the Mis, and in particular for the scalar 
amplitude Xi (d; s). 


4.6.4 The general case of four different masses 

Last but not least, we move to considering the most general configuration with four different masses. In 
this case the complexity increases even further and solving the IBPs in d dimensions brings to a reduction 
in terms of 11 different Mis 


Xi(d; s) = h{d-, 1,1,1,1,0, 0,0, 0,0), 
X 3 (d; s) = / 4 (d; 1 , 2 , 1 , 1 , 0 , 0 , 0 , 0 , 0 ), 
X 5 (d; s) = / 4 (d; 1 , 1 , 1 , 2 , 0 , 0 , 0 , 0 , 0 ), 
X 7 (d;s) = d 4 (d; 2 , 2 ,l,l, 0 , 0 , 0 , 0 , 0 ), 
Xg(d;s) = / 4 (d; 2,1,1,2,0,0,0,0,0), 
Xii(d;s) = h[ 


X 2 (d;s) =d 4 (d; 2 ,l,l,l, 0 , 0 , 0 , 0 , 0 ), 
X 4 (d; s) = / 4 (d;l,l, 2 ,l, 0 , 0 , 0 , 0 , 0 ), 
Is{d-, s) = hid; 3,1,1,1, 0,0, 0,0, 0), 
X 8 (d;s) = / 4 (d; 2 ,l, 2 ,l, 0 , 0 , 0 , 0 , 0 ), 
Xio(d;s) =/ 4 (d;l, 2 , 2 ,l, 0 , 0 , 0 , 0 , 0 ), 
:; 1 , 2 , 1 , 2 , 0 , 0 , 0 , 0 , 0 ). 


(4.64) 


The number of independent master integrals is obviously very large and, if all differential equations for the 
11 Mis were to be coupled, this would imply an 11-th order differential equation for any of the masters 
and in particular for the scalar amplitude Xi(d;s). It is therefore very interesting in this case to know 
how many Mis can be decoupled using the methods described above. Again it is enough to repeat the 
reduction to Mis, but fixing this time d = 2, and we immediately find that 5 out of the 11 Mis become 
linearly dependent and can be expressed in terms of the other 6 . Which Mis survive depends of course on 
the internal algorithm for the solution of the IBPs. In our case we find as independent Mis 


Xi(2; s) = hi2; 1, 1,1,1, 0, 0,0, 0,0), X2(2; s) = hi2; 2,1,1,1, 0, 0,0, 0,0), 

X 3 ( 2 ; s) = hi2; 1, 2, 1,1,0, 0,0, 0,0) , X 4 ( 2 ; s) = hi2; 1, 1, 2,1, 0, 0,0, 0,0), 

X5(2; s) = hi2; 1, 1,1,2, 0,0,0,0,0), Iei2; s) = hi2; 3, 1,1,1, 0,0,0,0,0) . 


(4.65) 


This implies that, the new basis of 11 d dimensional Mis defined following the recipe given above, fulfils a 
system of 11 differential equations, 5 of which decouple from the system as d —>■ 2. In this way we expect 
that a sixth-order differential equation for the scalar amplitude can be derived in d = 2 . 
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4.7 Comments and open questions 

Before moving on to the conclusions we would like to bring attention to some issues that might have gone 
unnoticed and which nevertheless leave room to very interesting open questions. In the previous sections 
we have worked out different applications of the ideas outlined in Section |31 We have seen explicitly that 
studying the IBPs in d = 2 or d = 4 can provide, for different Feynman graphs, identities useful to decouple 
the system of differential equations they fulfil. However in this discussion there is a point that we have 
avoided mentioning on purpose. Let us imagine to have to deal with a Feynman graph with three master 
integrals Ii, I 2 and I 3 , which fulfil a system of three coupled differential equations in the limit d —>■ 4, 
and let us suppose to apply the methods described above in order to try and decouple the system. We can 
imagine that by solving the IBPs in d = 4 only one relation can be found. In this case we know that we 
can use it in order to decouple one of the three integrals in the limit d —>■ 4, leaving therefore a system of 
two coupled equations, equivalent to a second order differential equation. Is this enough to say that there 
must be no other way to fully decouple all three differential equations? The answer is, in general, of course 
no. We have discussed already how the decoupling of the differential equations in any even integer number 
of dimensions is equivalent to the decoupling of the latter in d = 4. In Appendix |B] we show explicitly 
how, if one can find a basis that decouples the differential equations in d = 2 n, a corresponding basis can 
be constructed which decouples them in d = 4. Let us then go back to the problem of the three coupled 
master integrals. Let us imagine that studying the IBPs in d = 2 two linearly independent relations are 
found among the three masters and that this allows to completely decouple the system in d = 2. In this 
case we know that a corresponding basis would have to exist in d = 4 as well, and we could find it with 
the methods described in Appendix |B] On the other hand, if we had not found any new relations in d = 2, 
we could still decide to try in d = 6 , or in d = 8,10,12 etc. apparently without an end. Who or what tells 
us when to stop and, therefore, when we can assert without any doubts that not enough relations can be 
found, in any even number of dimension, to decouple completely the system? This question is extremely 
interesting and we unfortunately do not have a conclusive answer to it. For sure in all examples considered 
so far it has always been enough to study the differential equations in d = 2 and d = 4 only in order to 
find all needed relations to decouple the system. In those cases where not all equations could be decoupled 
(see for example Sections 14.31 14.41 14. 6 |) an attempt to consider different numbers of dimensions would 
simply produce no new relations at all, suggesting that there is no way to further simplify the problem, 
at least in this framework. Of course this does not constitute a mathematical proof in any respect. If 
these considerations have not brought to a definite answer yet, they nevertheless open the possibility for a 
different perspective in the way a system of differential equations for master integrals should be studied. 
We usually tend to think that the only physically relevant results are obtained when studying the system 
in the limit d —>■ 4. A lot of useful information, though, can be extracted studying the system as d —>■ 2 n 
and, sometimes, the relations found in this way appear to be independent and, in a sense, complementary. 
Whether these relations are really independent and how to determine the maximum amount of information 
that can be extracted by studying the IBPs in fixed integer numbers of dimensions remain open questions 
for now. It seems however reasonable to think that a more global approach, which allowed to study the 
systems of differential equations in general for any even number of dimensions (and not only in the limit 
d —> 4) could possibly bring a deeper insight in the structure and properties of the latter. 

5 Conclusions 

The method of differential equations has proven to be one of the most effective and promising tools for 
the evaluation of multi-loop and multi-scale Feynman integrals. The usual procedure consists in reducing 
all Feynman integrals to a basis of master integrals through integration by parts identities, then derive 
differential equations satisfied by the master integrals and finally try and solve them as Laurent expansion 
in (d — 4). For many problems of physical interests the coefficients of the Laurent expansion of the master 
integrals can be expressed in terms of a particular class of special functions called multiple poly logarithms. 
It has been noted that, whenever this is possible, a basis of master integrals can be found such that their 
differential equations become triangular in the limit d —>■ 4, allowing a simple integration of the differential 
equations by quadrature. It was moreover conjectured that in all such cases a canonical basis can be 
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found, turning the integration of the differential equations into a straightforward algebraic problem. If a 
complete set of boundary conditions is also known, the problem can therefore be considered as completely 
solved. On the other hand, different cases are known where such a simplification cannot be achieved and a 
minimum number of differential equations remain coupled. Of course in all these cases also a canonical basis 
(in the original sense introduced in [52]) cannot be found. Whenever this happens, it becomes of crucial 
importance to be able to determine the minimum number of master integrals which cannot be decoupled 
from the system. If two or more equations are coupled, in fact, no general technique exists to find a 
solution and one must resort to different considerations in order to find a complete set of homogeneous 
solutions, which can then be used in order to build up the inhomogeneous solution using Euler’s method 
of the variation of constants (see for example [39]). Of course, the larger the number of coupled equations 
is, the more difficult it becomes finding a complete set of solutions. Reducing the order of the system of 
differential equations is therefore essential from a practical point of view in order to be able to successfully 
tackle the problem. The issue is nevertheless interesting also from a more general point of view. Master 
integrals satisfying higher order differential equations, in fact, usually cannot be expressed in terms of 
multiple polylogarithms only and a very intensive theoretical effort has been recently devoted to determine 
the properties of the new special functions required. The most famous example is the two-loop massive 
sunrise graph. In this case two differential equations remain coupled and therefore the problem amounts to 
solving a second-order differential equation. It has been recently shown that the solution of the latter can 
be expressed in terms of a new generalisation of the multiple polylogarithms, called elliptic poly logarithms. 
Many questions are nevertheless still to be answered. Are elliptic polylogarithms enough for describing 
all Feynman integrals whose evaluation can be reduced to a second order differential equation? And what 
about higher order equations? 

A first step towards an answer to these questions seems therefore to be in a criterion to determine, 
given a set of master integrals and the system of differential equations they fulfil, the minimum number 
of differential equations coupled, and therefore the class of special functions required. In this paper we 
presented a simple idea which proved to be very useful in this respect. We showed in particular that the 
study of the IBPs for fixed integer values of the space-time dimensions, d = n, can provide the information 
required for decoupling the differential equations in the limit d ^ n. Indeed our criterion is, in principle, 
a sufficient but not a necessary one, in the sense that we did not prove that if no extra relation can be 
found among the Mis when d = n, then no decoupling is possible for d —^ n. The criterion has moreover 
proven to be extremely effective, inasmuch as it provided, in all cases that we considered, relations useful 
for decoupling some of the differential equations and therefore substantially simplify the problem at hand. 
It would indeed be extremely interesting to prove whether this criterion is not only a sufficient but also 
a necessary criterion, checking, for example, whether the number of independent Mis of the three-loop 
banana graph in d = 2 (see Section EM can be further reduced by any other means, reducing in this way 
also the maximum degree of the differential equation satisfied by the scalar amplitude. The criterion is 
moreover extremely simple to apply, since it can be very easily implemented into any existing public or 
private IBPs reduction code. In this respect, the possibility of pairing the study of IBPs in fixed numbers of 
space-time dimensions together with the new concept of a (pseudo-)finite basis of Mis, recently introduced 
in m, looks particularly promising. The latter, in fact, could potentially provide a way to automatically 
determine the highest poles developed by different Feynman integrals for different values of the space-time 
dimensions. 
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A Comparison with the Schouten Identities 

In this appendix we would like to show explicitly how the methods described in this paper are equivalent 
to the Schouten Identities introduced in |41j . We will consider again the two-loop sunrise with one massive 
and two massless propagators (see Section I4d1 for the definitions of the Mis) and try to derive relation (14.161) 
using the Schouten Identities. The two-loop sunrise graph depends on three independent momenta, the 
two loop momenta fc, I and the external momentum p. With these three momenta we can build up a 
d-dimensional Schouten polynomial which becomes zero as d —>■ n with n G N and n <2. Following [41] we 
start off by considering the quantity 

e(fc, l,p) = k^^rpP (A.l) 

defined in d = 3 space-time dimensions. Indeed (lA.ll) is nothing but the Gram determinant of the three 
vectors k, l,p. By squaring (IA.1|) in d = 3 we obtain a polynomial 

P2{d;k,l,p) = {ep,pk>^rpPf 

= P — k'^ {I ■ pY —p'^{k- lY — P {k ■ pY + 2{k ■ 1) {k ■ p) {I ■ p) . (A.2) 

The polynomial was obtained in d = 3 dimensions, but since it contains only scalar products of the three 
momenta it can be easily analytically continued to d dimensions and regarded as a d dimensional polynomial. 
By construction the polynomial is zero as d ^ 2 


P2(d ^ 2; fc,/,p) —>• 0 , i.e. P 2 {d ^ 2-, k,l,p) = 0{d — 2). 


Let us consider now the two master integrals defined in (14.31) . As discussed already both masters develop 
a double pole in (d — 2) 


Xj{d]p^) = ^ ^^(2;p2) + .^_r( ^)(2;p2 )+xf^( 2;p2) + , withj = l,2. 


(d-2)2 J ^ ^ ' (d-2) ^ 

We consider now the following quantities 


Z(d; 711,712, ns) = J Ti'^kTi'^l 


P2{d;k,l,p) 


(fc2)"i (^2)"2 + p)2 _ > 


(A.3) 


(A.4) 


which are of course linear combinations of integrals belonging to the sunrise graph, equation (14.11) . The 
Schouten polynomial goes to zero as d —> 2 and provides therefore an additional factor (d — 2) in the 
numerator, which can be used in order to alleviate the total divergence of the integral. Assume now that, 
using this piece of information, we can prove that the integral Z(d; ni, 712,713), for a given choice of the 
indices {711,712,713}, can develop at most a single pole in (d — 2) 


Z{d —>■ 2; Til, 712,713) c)c O 


1 

d^ 



(A.5) 


If this is true then, for any values of the indices {tii, 712,713}, we can consider the integral Z(d; ni, 712,713) 
as an integral of the sunrise family and use d-dimensional IBPs to reduce it to the two Mis 

Z(d; 711,712, 713) = (d; p2) I, (d; p2) + ^2) ^2) ^ g) 

where the are in general rational functions ofp^ and d. Suppose now that the C'"^"^"^(d;p^) 

do not develop any overall factor (d — 2) in the numerator. If this is the case, upon expanding in Laurent 
series both the left- and right-hand side of (IA.6I) and using the fact that the l.h.s. has only a single pole, 
we can find a relation between the double poles of the two Mis Ii and 12- To see how this works in 
practice, let us study this example for some specific choices of the indices 711,712,713 and see how this piece 
of information can be easily read out. 
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a) We start considering the easiest case ni = n 2 = = 1, i.e. we study the integral 


Z(d; 1,1,1) 




P2{d-,k,l,p) 

k"^ P {{k — I +p)^ — m^) 


Since /(d; 1,1,1,0, 0) has a maximum pole l/(d — 2)^, we could naively expect that, due to the overall 
(d — 2) factor carried by the Schouten polynomials, Z(d; 1,1,1) should develop at most a single 
pole l/(d — 2). This is in general of course not granted, since the polynomial in the numerator can 
worsen the UV behaviour of the integral. This naive expectation can be nevertheless easily verified by 
different means, for example using sector decomposition [^. On the other hand, as already specified 
above, the two Mis have both a double pole l/(d — 2)^. Expressing Z{d; 1,1,1) in terms of Ii and 
I 2 one finds easily that the corresponding coefficients C'i^^(d;p^), C' 2 ^^(d;p^) do not have any overall 
(d — 2) factor. Therefore we can expand both the left- and the right-hand side in Laurent series in 
(d — 2) and keeping only the first orders we get 


O 


d- 2 


(rrPY 


{d-2Y 


{ l [ ^\ 2 ; p ^)-{ p ^- m ^)4 ^\ 2 ; p ^))+0 


d-2 


(A.7) 


Eq. (IA.7I) gives for consistency a relation between the double poles of the two Mis, i.e. 


l[ ^\2;p^)-{p^-m^)4 ')( 2 ;/) = 0 , 


which is, as expected, identical to the relation obtained studying the IBPs in d = 2, Eq. (14.161) . 


b) As a second example, let us consider the case ni = 713 = 2 , 712 = 1 (or, equivalently, 712 = na = 2, 

Til = 1 ). Also in this case /(d; 2,1, 2,0,0) has a double pole in (d — 2) and we would naively expect 
that Z(d;2,l,2) should therefore develop again at most a single pole. This naive expectation can 
be once more verified explicitly using, for example, sector decomposition. Expressing Z(d;2,1,2) as 
linear combination of Mis and expanding in (d — 2) we get 


O 


1 


1 




^d-2j (d-2)2 

which for consistency implies 

_ (p2 _ ^2)j(-2)(2.p2) ^0^ 


d-2 


(A.8) 


in agreement with the previous case. 


c) As last example we can check what happens for the combination of indices tii = 772 = 2, 713 = 1. 

Again, repeating all considerations above, one finds that /(d; 2, 2,1, 0,0) has a double pole l/(d — 2)^ 
and Z(d;2,2,l) has instead only a single pole l/(d — 2). Reducing Z(d;2,2,l) and expanding in 
(d — 2) we find 

° (l^) = (d^ 't (^7*17") - (7“ - + O , (A.9) 

which once more implies 

l[-^\2-,p^) - {p^-m^)4-"\2-,p^) = 0, 
again in agreement with our previous hndings. 


Using the Schouten Identities one recovers therefore the same relation found by studying the IBPs in 
the limit d = 2. The relation can then be used, as shown in Section [4.1.11 in order to decouple the system 
of differential equations for this two-loop sunrise graph. We want to stress here that, differently from 
the direct study of the IBPs in d = n, Schouten identities can be derived only if a sufficient number of 
independent momenta exist. In the case of the two-loop sunrise studied here there are three independent 
momenta and we can therefore derive a Schouten identity in d = 2 only, but we have no mean to study 
possible relations among the Mis in d = 4. On the other hand, by studying the IBPs in fixed number of 
dimensions, one can easily try and look for relations in any number of dimensions, and in particular in 
d = 4, see Section l4T] 
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B Dimensional shift of systems of differential equations 

In this paper we showed that by studying the IBPs for fixed numbers of space-time dimensions d = n, 
one can in general find relations which allow to decouple some of the master integrals from the system of 
differential equations, and simplifying therefore the solution of the latter. In physical applications we are 
of course interested in the case c? = 4. On the other hand we have already briefly discussed how, once 
the full set of Mis is known as Laurent expansion for d -> 2 n, with n e N, by using Tarasov-Lee shift 
identities |431I44) one can reconstruct their Laurent expansion in any other even number of dimensions, 
and in particular in d = 4. It is therefore clear that, if by any means one can find a basis of Mis whose 
differential equations are in a convenient form as d —>■ 2n (triangular form, canonical form, etc.), then 
there must exist a corresponding basis of Mis whose differential equations look identical under the formal 
substitution (d — 2 n) (d — 4). In this appendix we want to show how this is indeed true and that such 
a basis can be obtained straightforwardly by a repeated use of Tarasov-Lee identities. 

Let us start considering a topology with N Mis Mj{d; Xij) where j = and we made explicit the 

dependence on the dimensions d and on the invariants of the problem Xij = pi ■ pj. Let us assume that, 
similarly to the case of the two-loop massive sunrise, Section 14.31 we are able to find a basis of Mis such 
that the differential equations take a particularly convenient form in the limit d —^ 2. In particular, in order 
to simplify the notation, let us assume that the differential equations are linear in d and can be written as 


■7^Mj{d;Xij) = Ao(2;xij) Mj{d;Xij) + {d - 2) Ai{2;xij) Mj{d;Xij ), (B.l) 

where Aq{ 2] Xij) is an N x N triangular matrix, which does not depend on d, while ^i(2; Xij) does not need 
to be triangular. If the system of differential equations is in this form, then by solving the homogeneous 
system (whose solution is now easier since the matrix is triangular) 


d 

— Hj (d, Xij ) — ^0 (2, Xij ) Hj (d, Xij ), 
oXfi 

we can find a transformation that puts the system of equations in the form 

d 


dx, 


-M'Jd; Xij) = (d - 2) B{x^j)MUd; Xij ), 


(B.2) 


(B.3) 


where Bn{xij) is an N x N matrix whose entries do not depend on d. Note that this does not ensure per 
se that the system will be in canonical form, since the entries of the matrix Bn{xij) might not be in d-log 
form, but could contain more complicated functions of the external invariants. This form is nevertheless 
very convenient for the explicit integration of the equations as Laurent series in (d — 2). 

In physical applications we are usually interested in the master integrals expanded in (d — 4), which are 
the physical space-time dimensions. Indeed we might think of solving the system (IB.31) as Laurent expansion 
in (d — 2), and then transport back the results to (d — 4) using Tarasov-Lee identities. Nevertheless we 
might also try and proceed differently and use Tarasov-Lee identities directly at the level of the differential 
equations in order to determine a new basis of Mis whose differential equations are identical to (IB.31) with 
the formal replacement (d—2) —?► (d —4). Tarasov shifting relations indeed contain this piece of information. 
By applying Tarasov shifting operators on the N Mis Mj (d, Xij) and reducing the result to the same set of 
Mis, we find N relations for the N masters which read 


Mj{d - 2, x^j) = ^ C[^\d\ x^j) Ml {d; Xij) 
i 




(B.4) 


where the cl^\d,kx) are rational functions of the dimensions d and of the invariants Xij. Define now a 
new set of N Mis Ij{d] Xij) as 

Ij{d-Xij)='^cj;^\d;x,j)Ml{d;Xij), j = l,...,iV. (B.5) 

i 
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Because of (|B.3I) . sending d ^ d —2, we find that 


-^—M'^{d-2;xij) =Ml{d-,Xij)^ = {d - 4) B^ix^j) M-{d - 2, Xij), (B.6) 

which can rephrased in terms of the new Mis as 

d 

-^^Xj{d]Xij) = [d - 4) B{xij)Xj[d\Xij). (B.7) 

Eq. (|R7| is exactly what we were looking for, namely a system formally identical to (IB.31) under the 
replacement {d — 2) —>■ {d — 4). While we showed this for a set of differential equations in a very special 
form (IB.31) . the considerations explained above are of course valid for any system of differential equations. 
Given a set of masters M = (Mi,..., M^r) and their system of differential equations in matrix form 

9 - 

-^^M{d;xij) = A{d;xij) M{d;Xij ), (B.8) 

where no constraint is applied on the matrix A{d] Xij), by using Tarasov relations we can define a new basis 

X{d;xij) = M{d-2-,Xij), (B.9) 

which, by construction, fulfils a system of differential equations in the form 

-^—X{d]Xij) = A{d- 2;x^j)i{d-Xij). (B.IO) 

OXij 

Of course, using Lee identities (or inverting Tarasov identities above) one can work in the opposite 
direction, shifting d ^ d4-2. Defining the new basis of Mis as 


J{d;Xij) = M{d + 2,Xij) , j = l,...,N, (B.ll) 

and following the same argument we find immediately 

d - 

-— J(d;Xij) = A(d + 2; Xij) J{d;Xij) . (B.I2) 

uXn 
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